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Abstract 


In May 1995, a set of spectrophotometric curves of the system Pluto-Charon was recorded with 
the UKXRT telescope equipped with the spectrometer CGS4. As for the previous observations by 
Owen et aL (1993), the spectra cover a part of the near infrared range, between 1.4 and 2.55 pm, but 
with a higher resolution of approximately 700. In both the 1992 (Owen et aL 1993) and 1995 data, the 
existence of solid methane is confirmed by numerous absorption bands, and the carbon monoxide and 
the nitrogen ices are identified by their respective signatures at 2.35 and 2.15 pm. The solid nitrogen 
seems to be the principal icy component and forms a matrix in which the CHU and CO molecules are 
diluted. However a spectroscopic analysis of the 1995 observations indicates that pure methane may 
coexist with its diluted phase in N 2 . In order to derive the horizontal and vertical distribution of these 
different species and to obtain some quantitative information about their characteristics, we have 
modeled the spectrum of May 15 that corresponds to the maximum of Pluto’s visible light curve. This 
was achieved by means of a radiative transfer algorithm dealing with compact and stratified media 
(Doutl and Schmitt 1998). Among the various representations we have tested to describe the surface of 
Pluto, only a geographical mixture of three distinct units explains all the significant structures of the 
analyzed spectrum. The first unit is a thin granular layer of pure CH 4 covering a compact 
polycrystalline substratum of N 2 -CH 4 -CO , which are in a molecular mixture (concentrations of CFL, 
and CO of the order of 0.45%, 0 . 1 - 0 . 2 % respectively). It covers about 70% of the observed area and 
corresponds to volatile deposits that are sublimating under solar illumination. The second unit is either 
(a) a single thick layer of pure granular methane Or (b) a unit similar to the first unit but with the two 
components inverted (i.e. with CH 4 forming a substratum and the N 2 -CH 4 -CO mixture a superficial 
layer of fine grains). Covering 20 % of the surface, it represents some old surfaces that have been 
sublimated for a long time, and eventually recovered later by very small amounts of fresh deposits of the 
molecular mixture N 2 -CH 4 -CO. Finally, the third unit may result from the condensation of very fine 
grains of nearly pure N 2 . It covers the remainder of the surface (about 10 %). All these results allow a 
better understanding of the processes of deposition, metamorphism, sublimation and transport affecting 
the different ices delected on PJulo during its climatic cycles. 
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1. Introduction 


Pluto, the ninth and the last identified planet of the solar system remains the only one that has 
not yet been visited by a spacecraft. However, our knowledge about this object has significantly 
progressed during the last few years thanks to recent instrumental developments. 

For the first time, the treatment of HST (Hubble Space Telescope) images, produced in late June 
and early July 1994 by the FOC (Faint Object Camera, ESA) in blue light has allowed Stem et al. 
(1997) to generate maps nearly covering the two hemispheres of Pluto with 160 km resolution. 
These maps show a huge spatial variability of the visible albedo with a dozen distinctive features, 
including a bright polar cap. The features are thought to be mainly related to changes of the 
chemical and physical properties over the surface. 

In the same way, the improvement of ground-based near-infrared (NIR) and to a less extent 
visible (Grundy and Fink 1996) spectrophotometry, has dramatically increased our capabilities to 
determine surface properties. This new generation of observations has come essentially from 
UK3RT (United Kingdom InfraRed Telescope, Mauna Kea, Hawaii). The initial measurements 
were made in 1992 with the NIR spectrometer CGS4 (Cooled Grating Spectrometer) and covered 
the 1.4 to 2.5 pm range at a spectral resolution of 6.7 nm (Owen et aL 1993). They included the 
light reflected by both the spatially unresolved disks of Pluto and Charon. A spectroscopic analysis 
of the different reflectance spectra obtained has provided crucial information about the nature and 
the physical state of the materials that cover the uppermost millimeters or centimeters of the surface 
of Pluto. The presence of solid methane was confirmed by eight characteristic absorption bands. In 
addition, carbon monoxide and nitrogen ices were identified by their respective signatures at 2.35 
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and 2.15 pm. Carbon dioxide, an ice that is detected on Triton, seems to be absent on Pluto, as well 
as hydrocarbons like C 2 H 2 , QHL,, C 2 H 6 , ... but water ice cannot be excluded (Owen et al. 1993). 

This study also showed that the methane molecules are most probably diluted in nitrogen ice 
but not as diluted as on Triton. This was deduced by noting that the position and width of most of 
the methane bands fit those of a laboratory sample of 1% CH 4 mixed in N 2 better than those of 
pure CH 4 ice on one side, or those of CR, even more highly diluted in N 2 , on the other side 
(Schmitt et al. 1992, Quirico 1995). The higher mean CH4 concentration (2.6 %) deduced by the 
model for Pluto compared to Triton also pointed toward an intermediate state of mixing. 
However, the laboratory data available at that time were insufficient to define this state more 
precisely. Subsequent extensive laboratory studies on pure CH4 and CO (as well as several other 
molecules) and their various mixtures with solid nitrogen (Quirico and Schmitt 1997a,b), Quirico 
et al. 1996, 1998, Schmitt et al. 1998) allowed a detailed spectroscopic analysis of the initial 

1992 UKIRT observation as well as a set of spectra taken in 1993. This led to the suggestion that 

tf 

a segregation of methane with an enriched top layer possibly occurs at the surface of Pluto 
(Schmitt et al. 1994, Quirico 1995) rather than a horizontal distribution of patches of pure and 
diluted CH 4 (Cruikshank etal. 1997, Grundy and Fink 1996). 

The analysis of the profile of the 2.15 |im band of nitrogen, used as a « thermometer » in the 

1993 observations, led Tryka et al. (1994) to propose a value of 40 ± 2 K for the temperature of 
nitrogen ice on Pluto. 

For Charon, water ice is the only compound that was recognized unambiguously by its 
spectral signatures at 1.55 and 2.04 pm (Buie et al. 1987), but other components may also exist 
(e.g. Roush 1994). 
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In May 1 995, an additional series of three Pluto-Charon observations was acquired with the 
UKIRT-CGS4 system in the same spectral range but with a spectral resolution improved by a factor 
of about 2.5 (2.64 nm). In this paper, we present a spectroscopic analysis of these observations as 
well as numerical modeling of one of these spectra. The scope of this work is a first attempt to 
derive the horizontal and vertical distribution modes of the different ices detected, in order to solve 
the problem raised above. This study also leads to more quantitative information about the intrinsic 
chemical and physical properties (composition, grain sizes, ...) and about the geographical extend of 
the solids identified so far. In the first part of the paper, the observational data is described as well 
as the relevant calibration procedures. The second part is dedicated to the analysis of the Pluto 
spectrum. First, we recall the contributions of the modeling approach, and present the framework 
and the numerical tool used for our simulations. Second, after modeling and extracting the 
contribution of Charon to the measured spectral reflectance, we assign the different bands observed 
and analyze the methane bands in terms of the physical state in which the methane occurs. We then 
describe how we tested and selected the best surface representations corresponding to the different 
coexistence modes of the molecular mixture N2-CH4-CO and pure methane suggested by the 
spectroscopic analysis. The synthesis of the previous results, followed by their interpretation and 
discussion is presented in the final part of the paper. 

2. Presentation and calibration of the observations 
In May 1 995, three non-spatially resolved reflectance spectra of the Pluto-Charon system were 
acquired respectively on May 11,13 and 15. The longitudes of the sub-earth point on the surface of 
Pluto were respectively 66°, 313° and 200°. 
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The CGS4 spectrometer measures the flux received by the UKIRT telescope with a spectral 
resolution of 2.64 nm (resolving power A/dX of approximately 750). The 1. 4-2.1 pm H range and 
the 1 .9-2.55 pm K range were recorded separately. The measurement of the lines of an argon lamp 
and the observation of the atmospheric bands of the OH radical, both in emission, ensure a good 
spectra] wavelength calibration with a mean uncertainty of ±1 nm. The details of the calibration 
procedure are described in Quirico et al. (1998). 

The radiative power which arrives on the detectors is integrated for one hour on the average. The 
resulting curve S(X), measured in digital units, is deeply marked by absorptions due to the transfer 
through the terrestrial atmosphere. As a consequence, a spectrum S«(X) of a solar type star is 
recorded in the same region of the sky with only a short time shift. Because S«(A) is affected by the 
same air mass, it allows us to correct these disturbing atmospheric bands (principally those of CO 2 
and H 2 0). Unfortunately, it does not insure a satisfactory radiometric calibration. For that, another 
reference measurement S ca ](A) is performed on a star whose near-infrared flux F cal (A) is well known 
but also disturbed by telluric absorptions. By the elimination from Scai(X) of the disturbing 
signatures which can be easily extracted in that case, one obtains its envelope S E cai(X) which is used 
for the calibration. The first quantity to be calibrated that way is the flux of the solar type star which 
can be written : 

Then dividing the radiance Bt(A) of a black body with a similar temperature by this flux determines 
an atmospheric correcting factor : 
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The equivalence between the solar type star and the corresponding reference black body can be 
written : 
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S*(A) being the envelope of the spectrum of the solar type star. The calibration operation 
continues with the multiplication of the spectral observed signal S(A) coming from Pluto-Charon by 
the correcting factor k(X), along with its transformation into a flux : 


pr.n- sw S - (A) F a) 
( ) s,wsL(A) F " ,( 


(4) 


Finally, the astronomical spectrum is converted into a geometrical albedo A g according to the 


expression 


, D 2 d 2 F(A) 
A * ( )_ R 2 F S (A) 


(5) 


where D and d respectively correspond to the heliocentric and the geocentric distance of Pluto- 
Charon, R their mean quadratic radius, F (A,) the corrected and calibrated flux as received on Earth 
and F S (A) the solar flux at 1 AU, as tabulated by Smith and Gottlieb (1974). 


After applying this calibration procedure, it appeared that no photometric coherence is achieved 
between the H and K spectral regions as judged from the overlapping 1. 9-2.1 pm range. In the K 
region all three measurements seem to give consistent albedo spectra while the H region gives 
unexpectedly high and variable values. On the May 11, clouds attenuated the signal from the 
calibration star (BS5384) while its H band spectrum was being measured. On May the 13 and 15, 
the skies were clear, but we observed the star at H through a partially opened dome because of 
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concerns that the starlight would saturate the array of the CGS4 spectrometer. Following the data 
reduction procedures described above, we multiplied the H band spectra from the three nights by 
scaling factors (0.75, 0.47, and 0.375 for the 11, 13 and 15 May, respectively) in order to match the 
H and K band spectra in their regions of overlap. We note that the scaled K band fluxes are 
consistent with those measured in earlier years when a wider slit was used and a better 
photometrical calibration achieved. For example, a K band spectrum measured on April 16, 1994, 
corresponds to the visible light curve minimum of Pluto at 96° longitude close to that of May 1 1, 
1995 (66°). On this curve, their corresponding visible albedos differ by a bit less than 5 % (Fig. 1), 
(Buie et al. 1997). It appears that the fluxes measured over the K region are very consistent and that 
the scaled flux of the H spectrum of May 1 1 , 1995, also fit most of the corresponding range of April 
16, 1994 quite well (Fig. 2). In addition, a H spectrum recorded on May 27, 1992 at 51° longitude, 
i.e. only 15° (and 4% visible albedo) away from the 05/1 1/95 observation, also fits the scaled 95 H 
spectrum within a few percents. Nevertheless, part of the 2 surrounding spectra does not match 
exactly with the 95 measurement, essentially outside the strong methane bands. The noticeable 
difference in observation conditions (-15° and +30° for the 92 and 94 spectra relative to the 95 one) 
as well as a possible real small evolution of the surface of Pluto in terms of average surface 
composition and/or physical state, may explain the spectral differences. 

In conclusion, we believe that our initial photometric calibration is correct for all three spectra in 
the K region and after rescaling for the H region spectra. Knowing that the 1992 and 1994 reference 
spectra at lower resolution display an estimated uncertainty of ±15%, we think that our calibration 
operation yields all three H&K May 1995 flux spectra within ±20%. In the following we will use 
the calibrated and concatenated 1.4-2.55 pm flux spectra, except for May 13 restricted to the H 
region due to some spurious effects in the K region data. 
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3. Analysis of the observation 


3.1 Framework of the modeling approach 

In the visible and near infrared ranges, the flux that comes to the earth from the Pluto-Charon 
system or from any other planetary body in general, totally or partly arises from the reflection and 
the absorption of the solar light which penetrates more or less deeply into the surface. The intimate 
interaction with the constituent materials implies that their physico-chemical and structural 
properties potentially mark the spectral signal. Moreover, as light attenuation into the materials 
often changes by several orders of magnitude at different wavelengths, every reflectance level of a 
spectrum probes a particular range of depths. If it is modulated enough, the general shape of the 
spectrum then becomes very sensitive to the occurrence of a stratification of the surface properties. 
For all these reasons, beyond a simple identification of the chemical nature and of the physical state 
of the components, as expressed by the precise absorption band positions (see part 3.3), the 
modeling of the intensity and the shape of the bands is the only efficient way to obtain additional 
information (see part 3.4). These are mainly the chemical composition of the compounds and the 
texture, as well as the granularity of the materials they form. These properties can be obtained as a 
function of depth in the case of surface stratification. 

However things can be a bit more complicated. Most of the time different kinds of uniform surface 
units, which form something like a more or less regular geographical mosaic that is not spatially 
resolved by the observation, may contribute to the same spectral signal. Then, they linearly add 
their respective contributions according to their surface proportion and the viewing geometry. This 
kind of addition is called a 'geographical mixture'. This is always the case for a spectrophotometric 
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measurement, like the one we analyzed, because it integrates the light reflected by the whole 
hemisphere of the planetary body of interest. 

For the modeling of the spectrum of the Pluto-Charon system, we used an efficient radiative 
transfer algorithm calculating the spectral bidirectional reflectance of an icy or mineral stratified 
surface (Doute and Schmitt 1998). Each layer of such a surface can present a granular or a compact 
texture. The aim of the algorithm is to determine, for any kind of illumination and viewing 
conditions, the link between the physical, chemical and structural properties of the underlying 
medium and its observed spectral signature. A potential stratification is described by the vertical 
superposition of any number of plane-parallel homogeneous layers. 

Practically, the algorithm is integrated into a software that can globally simulate and then analyze 
by interactive inversion hyperspectral data according to a wide variety of observational and surface 
situations (Doute 1998). Indeed this system can deal with spectrophotometric data (as is the case 
here with Pluto-Charon), and with hyperspectral images, both corresponding to surfaces that can 
present any kind of chemical composition, geographical and vertical distribution of the materials. 
The modeling is performed by the integration of the radiative spectral power reflected by the 
different surface units that generally constitute a planetary surface. The previous radiative transfer 
algorithm ensures the calculation of each contribution of such a geographical mixture. All the above 
mentioned characteristics lead to a very powerful system to reach our main objective : to test and 
select different surface representations of Pluto. 

3.2 Modeling and extraction of the contribution of Charon 

We first need to separate the respective contributions of Pluto and of its satellite. The problem 
consists of extracting a spectrum of Charon from our observation which was acquired at a far 
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higher resolution (dX = 2.64 nm) than any available Charon spectrum. Roush et al. (1996) derived 
the spectral geometric albedo of Charon from very poorly resolved (dX = 0.1 pm) 
spectrophotometric data acquired between 1.5 and 2.5 pm during a special mutual event : the total 
eclipse of Charon by Pluto on April 23, 1987 (Buie et al. 1987). We first modeled this geometric 
albedo curve with high resolution simulations deresolved at dX = 0.1 pm (Fig. 3). For that, we have 
to choose a suitable, but simplified, surface representation of Charon. The simplest one we found is 
a uniform hemisphere covered by a geographical mixture of two kinds of surface units : 

• a unit U w , which is a granular optically semi-infinite layer of water ice, 

• a unit U n , which is spectrally neutral and of unknown chemical nature. 

For the first unit, intended to reproduce the absorption bands, we use new optical constants of 
crystalline water ice at 40 K (Grundy and Schmitt 1998). The sedond unit permits the 2 pm flux of 
the strong H 2 0 bands to rise to the correct level. This was impossible to realize with a hemisphere 
entirely covered with H 2 0 grains. The principal parameters inferred by the inversion (surface 
proportion of the water ice ^“^20 = 60 %, mean grain size d H2 o = 50 pm) as well as the chosen 
surface representation may not be very significant, nor well constrained, but the deresolved 
spectrum they generate by direct modeling is within the error bars of the geometric albedo of 
Charon. 

Using the previous values, we calculate the corresponding high resolution (dX = 2.64 nm) spectrum 
of Charon which is directly transformed into flux using Eq. 4. This flux contribution is then 
removed from the Pluto-Charon 1995 observations. We finally convert the resulting spectrum of 
Pluto back to geometric albedo thanks again to Eq. 5 (Fig. 4). 

As the eclipse observation of Charon corresponds to a Pluto central longitude of 190°, the removal 
of its spectral contribution is strictly valid only for the May 15, 1995 observation (200° Pluto 
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longitude). Near-infrared spectral variations with Charon longitude have been observed and 
attributed to changes in water ice absorption (Bosh et al. 1992, Buie and Shriver 1994), so the 
removal of the eclipse observation from the May 11 or 13, 1995 observations may lead to incorrect 
subtraction of the water ice signature from Charon. For that reason, in this paper we restricted our 
modeling effort to the May 15 observation. On the other hand, the following spectroscopic analysis 
of the two other Pluto spectra (at 66° and 313°), mainly based on CH4 band positions, should not be 
significantly altered by the first order removal of Charon water ice absorption using its 190° 
synthetic spectrum. 

3.3 Spectroscopic analysis 

After some discussion of band assignments we will present the spectral analysis applied to 
the higher resolution spectra described in this paper. 

3.3.1 Band assignment 

The high resolution spectra of Pluto exhibit 20 bands in the 1. 4-2.5 pm range. About half of 
them were observed and assigned by Owen et al. (1993) to N2, CH4 and CO. We complete this 
identification and indicate, for each band, the specific vibration mode of the corresponding 
molecule (Fig. 5a,b). In particular, the observation of the second overtone of solid CO at 1.579 
pm confirms the identification of this molecule on Pluto. The band at 1.69 pm that we attribute 
to pure solid CH4 is a special case discussed in more detail below. 

Compared to Triton, there is no positive sign of the 3 strong CO2 ice bands around 2 pm, 
thus excluding the presence of this molecule in significant amount at the surface of Pluto. The 
problem of water ice is more complicated because its broad absorption bands combined with the 
broad CHL, bands make its direct identification very difficult. As an illustration, just looking at 
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the spectrum of the Pluto-Charon system displayed in Figure 4, it is not evident that it contains 
the signature of water ice, although this component actually contributes to at least 13% of the 
projected surface! This point will be further addressed by the detailed modeling of the spectrum 
(see 3.4), but only for the observation simulated in this paper (5/15/95). 

Only one band that is present in all 3 of the May 1995 observations, but not observable in the 
former lower resolution observations, cannot be attributed to one of the three molecules 
identified so far on Pluto. It peaks around 1.75 pm, and is slightly shifted from one observation 
to the other in the 1.749-1.752 pm range. This band has also been observed in the spectrum of 
Triton where it has been suggested that ethane (C 2 H6), or perhaps propane (CjHs) diluted in solid 
nitrogen may be potential candidates (Quirico et al. 1998). Comparisons with the available 
laboratory data on various molecules (Quirico and Schmitt 1997a, Quirico et al. 1998) show that 
pure C 2 H 2 , C 2 H4, C 2 H 6 , C 3 H 8 , HCN, HC3N, CH 3 OH, SO* NO, N0 2 and NH 3 ices cannot explain 
this band but that C 2 H 6 and possibly C 3 H 8 diluted in solid nitrogen display a strong band only 
shifted by 0.003-0.006 pm from the observed feature. However, the limited relevance of the 
currently available data on both molecules (measured at 1% concentration in a-N 2 at 2 IK) does 
not explain the remaining shift, thus still precluding a firm identification. Laboratory data on 
these molecules at more relevant temperatures and concentrations are clearly needed. 

3.3.2 Physical state of the CH 4 molecules 

Using the laboratory data of Quirico and Schmitt (1997a), we can compare the methane 
bands in the Pluto spectra with those of pure solid methane around 40 K and those of CH4 diluted 
in solid (3— N 2 around this same temperature (Fig.. 6a,b). The first striking point is that the weak, 
but real, band present in all 3 observations around 1.69 pm is only active in pure solid methane. 
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Indeed, this band has never been observed in any sample of CR, diluted in N 2 (neither in a or (3- 
N 2 phases) and at any temperature. This band is clearly a witness of the presence of pure methane 
at the surface of Pluto. 

The second point is that the peak of the strongest bands in the K region (at 2.20, 2.32, 2.37, 
and 2.43 pm), particularly for the 15 May spectrum, are situated at positions intermediate 
between pure and diluted methane, contrary to the same absorptions in the spectrum of Triton 
that perfectly fit the positions of diluted CH4 (Fig. 4b). These bands are also much wider, partly 
due to saturation of the bands (at 2.32 and 2.37 pm) but possibly also due to the simultaneous 
contribution of both forms of methane. If we now look at the medium bands in the H region, they 
more closely fit the wavelength of diluted methane for the 11 May observation, but they are still 
slightly shifted towards low frequency compared to Triton in the 15 May, and particularly in the 
13 May spectra (Fig. 6a). In addition the side band near 1.78 pm is less pronounced than on 
Triton, indicating a possible contribution of pure methane, partly filling the trough between the 
side and the main band. Finally the two weak bands around 1.48 and 1.95 pm follow similar 
trends, the 11 and 15 May bands being at intermediate positions and the 13 May bands almost 
fitting the positions of pure methane. 

All this information points toward the coexistence of pure and diluted methane at the surface of 
Pluto. For the 15 May observation the apparent increase of the diluted character of methane when 
looking at bands with decreasing strengths, probing deeper into the surface, also suggests that at 
least part of the pure methane should be preferentially accumulated near the surface. The reverse 
seems to occur for the 1 1 May spectrum, pure methane showing up at depth (weak bands) and 
diluted methane dominating the lop layers. The effective occurence of such vertical stratifications. 
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and of any other type of CH 4 segregation, should be confirmed by a detailed modeling of the whole 
spectrum for various possible surface representations (see 3.4). 

It should be noted that the phase diagram of the binary N2-CH4 mixture (Prokhvatilov and 
Yantsevich 1983) tells us that absolutely pure methane ice cannot coexist with nitrogen ice, but that 
a small fraction of nitrogen (3-4% around 40K) should be diluted in it. The positions of the CH 4 
bands in such a system are slightly shifted from the positions for pure CKL, ice, as shown by 
laboratory spectra (Quirico et al. 1996, Schmitt et al. 1998). However this shift is smaller than the 
resolution of the observations and does not affect our interpretation. In the following we will 
continue to refer to this Ohrrich mixture containing a few percent nitrogen as “pure methane”. 

These intermediate positions of the CH4 bands provide important information about the physical 
state of CH 4 at the surface of Pluto. On the other hand, they preclude the extraction of any 
information or constraint on the temperature of the N 2 -CH 4 -CO mixture and on the crystalline 
phase of solid N 2 based on the precise position of the CH 4 bands, as was possible for Tnton where 
only diluted CH 4 is present (Quirico et al. 1998). In particular, the most sensitive band relative to 
the phase of nitrogen is strongly perturbed by the suspected contribution of pure CPL, ice. The 
detailed modeling of the shape of the CH4 bands may help constrain the surface temperature, but it 
is anticipated that the blending between the two Oh, contributions will dominate the spectrum. The 
CO bands being weak and quite insensitive to temperature and phase, the only remaining 
temperature tracer seems to be the weak nitrogen band at 2.15 pm. As has been noted before, its 
broad shape is clearly indicative of the predominance of the P phase, i.e. a temperature larger than 
35.6 K. But only a thorough modeling of the shape of the weak N 2 band together with the 
underlying CH4 band may provide some additional constraint on the temperature of the nitrogen- 

rich deposits. 
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3.4 Test and selection of the best surface representations 


3.4. ] Presentation 

The previous spectroscopic analysis suggests the coexistence of two phases of methane on the 
surface of Pluto : the pure solid form and methane diluted in solid nitrogen. Three simple potential 
modes of coexistence can be imagined: (I.) a granular (intimate) mixture, (II.) a vertical 
stratification, and (HI.) a geographical mixture. An intimate mixture of grains respectively formed 
by pure methane and « CFL, contaminated » nitrogen is the first possibility (L). A second possibility 
(II.) would be a vertical segregation of pure CR» and of a N 2 -CH 4 -CO molecular mixture. The third 
possibility (in.) is that more or less extended horizontal patches of pure and diluted CH4 coexist on 
the surface of Pluto. Nevertheless, as shown in 3.3.2, these basic surface representations will have 
to be combined in different ways in order to provide the best solutions for the problem of the 
species distribution. Due to the geographically integrated nature of the available spectra, it is not yet 
possible to derive a precise spatial distribution of the different physico-chemical units we can 
expect at the surface of Pluto. Consequently, we propose a crude idealization of the three basic 
representations of differentiation we have described, each covering the entire visible hemisphere of 
Pluto : 

I. First basic representation Ri (intimate) 

A uniform deposit composed of one single thick layer : 

- a unit Ui formed by a homogeneous population of pure CH 4 and N2-CH4-CO grains 
intimately mixed (G c and G m , respectively), 
n. Second basic representation Rs ( stratified ) 

A stratified and uniform deposit LU, composed of two superimposed distinct layers : 
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- a superficial layer U of finite depth, formed by a homogeneous population of pure CH 4 
grain s. 

- an optically semi-infinite substratum L™ formed by a homogeneous population of NrChL,- 
CO (molecular mixture) grains. 

HI. Third basic representation Rg (geographic) 

A uniform geographical 'mosaic' of two distinct units : 

a unit U c formed by an homogeneous population of pure CH 4 grains. 

- a unit U m formed by an homogeneous population of N 2 -CH 4 -CO (molecular mixture) 

grains. 

These three configurations imply three different lists of free parameters for the models, which are 
listed in Table I. In the previous and all the following notations of the « Units », U, the subscript 
<<c » refers to pure CH4, « m » to the N2-CH4-CO molecular mixture, « m+c » (aliased to « i » for 
the « Representations », R) to the intimate (granular) mixtures of « c » and « m » grains, «c/m » 
(« s » for R) to the superposition of « c » over « m », « m/c » (« s » for R) to the reverse 
stratification, and « n » refers to the neutral or fine grained N 2 unit All surface representations 
(except the basic representation R g = Rem) that is a geographical combination of two or more of 
these units will be described by the concatenation of the corresponding subscripts (e.g. : Rjsn - Uc+m 
+ Uc/m Un)» 

The first difficulty that arises when modeling the spectrum of Pluto with either basic surface 
representations Ri, R s , or R g , is an obvious impossibility to fit the albedo level for the 2.32 and the 
2.37 pm CH 4 bands. Indeed, due to the required saturation of these two bands, the models always 
lead to albedos much lower than observed. There are two possibilities to increase the albedo. 
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- add a new geographical unit that will contribute most of the flux at 2.32 and 2.37 pm. This 
implies a minimum level of brightness for the constituent materials. 

- mix at least part of the surface grains with some other particles that have to be bright enough (high 
single scattering albedo) to increase the reflectance around 2.32 and 2.37 pm in particular. 

In both cases, the additional necessary component(s) should have an almost neutral spectral 
reflectance over the entire range, since we do not observe any significant continuum slope in the 
spectrum and only a single weak unidentified localized signature (possibly linked with units 
containing CH4, if the ethane identification is confirmed). Among the molecules that have been 
identified in the spectrum of Pluto, only nitrogen, even a with small amount of CH4 and CO, can 
present this neutrality, if it is present as grains that are sufficiently fine. This property comes from 
the weakness and the uniqueness of the N 2 absorption band at 2.15 pm. We choose to add a new 
unit of fine grained nitrogen ice to the geographical mosaic with the other compounds rather than 
forming a granular mixture of N 2 -rich grains within the unit U m , as it is thermodynamically very 
improbable to have an intimate mixture of two very different grain sizes of the same volatile 
material. In addition, because of the non-linear behavior of a granular mixture, it is necessary to 
introduce an extremely large quantity of fine grained nitrogen to produce a substantial rise of the 
reflectivity at 2.32 or 2.37 pm where methane is highly absorbing, whatever its state of dilution 
(pure or diluted). 

Consequendy, we add to all our future representations based on Rj, R s and R g a spectrally neutral 
geographical unit U n that is described by the following parameters : 

U n [fine grains of nitrogen] 

• surface proportion c P surf „ 

• mean grain size r n 
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scattering parameter g n . 


One can note that several scattering parameters ‘g’ appear in the previous list and in Table I. 
They express the angular distribution of the scattered power by a volume element of a layer and are 
related to the Henyey-Greenstein phase function as defined and illustrated by Hapke (1993). They 
Cannot really be constrained by a directionally integrated measurement as is the case for our 
spectrum of Pluto. Without any additional information, we decided to set all these parameters to 
zero (isotropic scattering grains). This particular choice may have a noticeable effect on the values 
of the parameters closely related to g, i.e. mainly the grain sizes, but not on their orders of 
magnitude. For the other parameters, the results of the inversion should be relatively insensitive to 
our choice. 

In order to calculate the scattering and the absorption properties of the materials that constitute 
the different layers of the surface units, the radiative transfer model we use for the analysis needs 
their optical indices versus .wavelength. The spectroscopic analysis has given some information 
about the physical state of the identified components : N 2 , CH4, and CO (see part 3.3). 
Consequently, we select a set of optical indices, measured in the laboratory and compatible with 
these constraints. Below, we describe these data, ice by ice : 

• Pure CH 4 ice : we use a set of measurements performed at 40K (unpublished) and similar to the 
data at 24 K and 30K (Quirico and Schmitt 1997a, Schmitt et al. 1998). The effects of 
temperature principally affect the position and the width of the bands (as well as the 
continuum), but remain weak in these data, especially when considering the resolution of the 
Pluto observations. Thus, the coherence of the optical index file is insured if the effective 
surface temperature of CR, ice is not loo different from 40K. This may be the case because 
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temperatures slightly higher (up to 54K) than the average are expected for this kind of deposit 
due to its lower volatility than N 2 ice (Stansberry et al. 1996). 

• Pure N 2 ice : the laboratory data we use were obtained in the P phase at 36.5K, slightly above 
the a-P phase transition, by Grundy et al. (1993) and Quirico 1995. 

• Molecular mixture NrCHLrCO : it is not possible to experimentally obtain the optical index 
spectra corresponding to all the potential values for the concentrations Cch 4 and Cco in solid 
nitrogen. As a consequence, we numerically generated the data needed from a sequence of 
experimental measurements concerning N 2 -CH 4 and N 2 -CO mixtures at 36.5 K obtained by 
Quirico and Schmitt (1997a,b). For that purpose, we considered three real or virtual 'pure' 
reference components. The first one is identified with pure N 2 ice, the second and the third ones 
with methane and carbon monoxide initially diluted in the nitrogen matrix, but artificially 
normalized to unit concentration (see Quirico et al. 1998). Then we assume that, for each 
wavelength, the scalar product of the complex indices of N 2 , of the normalized diluted CEL, and 
of the normalized diluted CO by the vector [(1 - Ccim * Cco). Cch4. Cco] gives the optical 
constants of the mixture. The normalization and scalar product operations do not change the 
spectral effects linked to the dilution. Indeed the slight high frequency shifts of the absorption 
bands of CHL, diluted in P N 2 compared to pure CEL, (4 to 6 nm in the H and K regions) strongly 
determine the modeling of the astronomical spectrum of Pluto. 

The inversion mode we use is an interactive trial-and-error sequence of direct modeling that 
progressively improves the fit of the spectrophotometric curve of Pluto. For that purpose we try to 
reproduce the main and more significant spectral characteristics of the spectrum. These spectral 
characteristics are the shape, the intensity and the position of the principal absorption bands of 
methane, nitrogen and carbon monoxide, as well as the presence or the absence of very distinctive 
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weak structures, all related to CHL, in its two phases. Aside from the bands of N 2 at 2.15 pm and CO 
at 2.35 pm, the methane features can be arranged for physical reasons into four independent groups 

(see Fig. 5a, b.) : 

. group one : the 4 strong bands of the K region at 2.2, 2.32, 2.37 and 2.42 pm. This group also 
contains the slight bulge around 2.28 pm standing on the high frequency wing of 2.32 pm band. 
This feature appears to be a significant structure because it can be reproduced by the modeling. 

• group two : the 3 medium bands of the H region at 1.67, 1.72 and 1.79 pm with, on the left 
wing of the first and last ones, a small side band with a position that is typical of methane 

diluted in nitrogen. 

• group three : the weak band at 1 .85 pm (H region). 

. group four : the little absorption of pure methane which appears at 1.69 pm on the hill 
separating the bands at 1 .67 and 1 .72 pm. 

3.4.2 Analysis 

The first task to achieve in order to precise the horizontal and vertical distribution modes of the 
different detected ices is to confirm one of the main results of the spectroscopic approach : there is a 
methane segregation at the surface of Pluto, whatever its nature. For that, we try to simulate the 
Pluto spectrum with even simpler representations than the three basic ones (R„ R. and R,), t.e„ 
considering only one form of CM,, either pure or diluted. 

Hie first R„, is a geographical mosaic of two units : one unit U„ formed by a molecular mixture 
of N 2 -CH,-CO and one unit U. (Nr). As is illustrated in Ftg..7, the simulations do not reproduce the 
positions of the methane bands that are intermediate between the pure and the diluted forms. 
Moreover it appears impossible to simultaneously fit, even roughly, the intensity and the width of 
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the N 2 and all the CH 4 bands. Finally the characteristic signature of pure CH 4 at 1.69 |im is absent. 
In a second simple representation, R^, the unit U m is replaced by an unit U c composed of pure CR» 
grains. Even more critical misfits are obtained (shifts, shapes, much too strong 1.69 pm band, ...) 
(Fig. 7). These models prove more directly that both pure and diluted methane coexist at the surface 
of Pluto, but they did not tell which types of differentiation occurs. 

In order to discriminate among the different basic types of differentiation, we performed the 
second task dedicated to test a series of more or less complex surface representations taking into 
account these differentiations. They are divided into three families. The first and the third families 
are respectively composed of the granular mixture R; and the geographical distribution R g , 
completed with the neutral unit U„ as stipulated before. The representations of the second family 
contain U n and at least one stratified unit (generally Uc/ m ). A given representation tries to reduce the 
limi tations of the previous one. In this paper, we choose to focus on the best fits of the Pluto 
spectrum we have obtained so far, with each representation. In each case, we will give a synthetic 
scheme of the representation, the list of the associated modeling parameters and their best fit values. 

a.) ‘Geographical representation R in with a granular mixture’ 

According to the phase diagram of the binary N2-CH4 molecular mixture (Prokhvatilov and 
Yantsevich 1983), the coexistence of almost pure CH4 grains (with small amount of nitrogen) and 
nitrogen-rich grains saturated by CH 4 can occur under some temperature and mean composition 
conditions. If the conditions are met, the concentration of CH4 in the N 2 matrix is constrained to a 
value of approximately 2% around 40 K. On the other hand the CO concentration is not constrained 
by the N 2 -CO phase diagram. 
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Taking these constraints into account, our simplest representation Rj n includes an unit U n and an 
intimate (granular) mixture of pure CHL, ice grains and N2-CH4-CO grains with the CH4 
concentration fixed to 2% (unit U c+m ). We took the optical constants of pure CH4 for the former 
grains, as only very small spectral changes are expected and measured for methane ice with 1-2% 
nitrogen diluted in it (Quirico et al. 1996; Schmitt et al. 1998). For the 2% CH4-N 2 (-CO) grains we 
used the values of diluted CH4 (stated limit of validity, see 3.4.1.). 

The best fit of the observation by this representation is shown in Figure 8. Although a 
satisfactory fit of most of the CH4 bands is achieved, there are at least three misfits that cannot be 
solved with this representation. First the strength of the N 2 band is more than an order of magnitude 
too weak. Increasing its strength by increasing its grain size strongly irremediably saturates the 
methane bands. The high (2%) CH4 concentration in nitrogen is surely the principal reason . 
Secondly, the low frequency wings of the strongest GHL» bands, especially at 2.32 and 2.37 pm, are 

not well reproduced, indicating some lack of pure methane close to the surface. On the other hand, 

< 

this representation has too much pure CH4 at some depth, as witnessed by the too strong 1.69 pm 
band. All these results, obtained with a vertically homogeneous distribution of the different types of 
grains, point to the existence of a stratification favoring the presence of pure methane close to the 
surface. This leads to the next type of surface representations. 

b.) ‘Geographical representations with stratification’ 

Among this category, the selected representations of the surface of Pluto were inspired by two 
preliminary trials. The first one R sn consists of a uniform geographical mosaic of two distinct units : 

• a unit Uc/ m associated with the basic vertical distribution R s , (a layer of pure CR» grains over a 
N2-CH4-CO mixture), 
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• a neutral unit U n composed with fine grains of nitrogen. 

The second trial, R s ' n , is similar to the first one but the two layers of the unit Uc/ m are inverted 
giving Um/c- Both simulations reproduce the main features of the spectrum quite well but do not 
show all the finest details at the same time (the bulge around 2.28 pm, the side bands in the H 
region and the absorption of pure methane at 1.69 pm). Nevertheless they tell us that, if a 
stratification of the deposits really occurs, pure methane should he at the surface but also at some 
depth. This conclusion leads to the first retained representation : 
i ) Representation R scn 

A uniform geographical mosaic of three distinct units : 

• a unit like Uc/ m of pure CH4 grains stratified over a N2-CH4-CO molecular mixture, 

• a unit U c formed by a homogeneous population of pure CH4 grains, 

• a unit U n of fine spectrally neutral grains of nitrogen. 

Table n contains a synthetic scheme of this representation, the list of the associated parameters as 
well as the values derived from the best modeling of the spectrum. Figure 9. illustrates the quality 
of the fit obtained between the observed and the synthetic spectra (see also Fig. 12a,b). 

In the K region (Fig. 12b), the spectral position, the width and the shape of the methane bands at 
2.20, 2.32 pm and to a lesser extent at 2.37 and 2.43 pm are quite weU reproduced. In particular, we 
can notice that the slight bulge around 2.28 pm is fit well by the synthetic spectrum. The absorption 
of nitrogen at 2.15 pm and of carbon monoxide at 2.35 pm correctly mark their presence by two 
small spikes. However, the asymmetry that affects the wings of CO band is reversed. If now, we 
move to the H region (Fig. 12a), the simulated spectrum reproduces the global shape of the 
intermediate bands of CH 4 , with the exception of the intensities of the 1.67 and 1.72 pm bands 
which are slightly too strong. The small side bands around 1.65 and 1.78 pm clearly stand out and 
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coincide rather well with the observation. On the other hand, we can notice that the small 
absorption spike of pure methane at 1.69 pm is present but a bit too weak. The shape and 
particularly the width and intensity of the 1.85 pm band are not well fitted, but between 1.82 and 
1.95 pm imperfect cancellation of telluric absorptions probably occurred in the spectrum of Pluto. 
Finally the continuum around 1.5, 2 and 2.5 pm is not well reproduced, although the 4 weak CFL, 
bands at about 1.45, 1.48, 1.93 and 1.97 pm are satisfactorily simulated. However we know that in 
these ranges of very weak absorption, our optical constants of CH 4 diluted in nitrogen are very 
uncertain and noisy (see Quirico et al. 1998). 

ii ) Representation R„~„ 

For the second selected representation, we consider a uniform geographical mosaic of three 
distinct units : 

• a unit Uc/ m (CtL* stratified over a N 2 -CH 4 -CO molecular mixture), 

• a second unit TU similar to the previous one but with the layers UandL™ inverted (N2-CH4- 
CO over pure CHLj), 

• a third unit U„ of fine spectrally neutral grains of nitrogen. 

Table in displays a synthetic scheme of the R ss ' n representation, the list of the model parameters as 
well as their best fit values. Figure 10. shows the best fit obtained between the observed and 
simulated spectra (see also Fig. 12a,b). 

There are very few differences between this simulation and Rsc„, whose qualities seem to be 
globally equivalent. We only notice, in the H region, the improvement of the intensity for the 
intermediate methane bands and the unfortunately correlated degradation of the intensities of the 
two side bands. The fit of the hill around 1 .75 pm also deteriorates a little bit 
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It should be mentioned that, among the different parameters which are used in the current 
representation, the concentrations (few and C 2 co as well as the mean grain size r^ of the molecular 
mixture superficial layer of the unit Um/ C are not well constrained. Indeed, this layer has a limited 
influence on the weak N2 band, which is essentially produced and controlled by the substratum of 
the unit Uc/ m - We can be much more confident of the parameter values obtained for the last- 
mentioned unit Uc/m- 

We will not proceed further with the stratified surface representations and keep the two best fitting 
simulations we have obtained so far. 

c.) 'Geographical representation R mcn without stratification’ 

In order to test if a stratification really has to be invoked, we first study a simple geographical 
distribution of pure methane and methane diluted in N2 (with CO). 

We thus consider a uniform geographical mosaic of three distinct units : 

• a unit U m formed by a homogeneous population of N2-CH4-CO (molecular mixture) grains, 

• a unit U c of volatile deposits made of a homogeneous population of pure CH4 grains, 

• a third unit U n of fine spectrally neutral grains of pure nitrogen. 

In fact, Rmcn correspond to the basic representation R g plus the neutral unit U„. The synthetic 
scheme of the representation, the list of the associated parameters as well as the best fit values are 
contained in Table IV. Fig 11 . illustrates the quality of the fit reached between the observed and 
synthetic spectra (see also Fig. 13 a,b). 

The global quality of this simulation comes close to the level of precision displayed by the 
previous modeling R scn and R ss ~ n but without reaching it completely. If we exclusively concentrate 
on the noticeable differences between R ss ~ n and R m cn, we can say that the latter fils the troughs on 
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each side of the 2.42 pm band of methane better but not the strong saturated absorptions at 2.32 and 
2.37 pn., which are too broad (Fig. 1 3b). In addition, the band at 2.32 pm does not coirecdy display 
the characteristic bulge around 2.28 pm, which is a significant feature. The signature at 2.2 pm is 
also a bit too strong and wide. If we now pay attention to the nitrogen and the carbon monoxide 
bands, we still give preference to the R„% simulation. Indeed, the IU simulation carves the former 
and lowers the latter a bit too much. On the other hand, the spectrum Rmen more longer fits the 
observation toward the continuum below 2.14 pm than the representation R„ In the H region 
(Fig. 13a), the agreement of the Rm™ simulation appears to be significantly less satisfactory if one 
refers to all the significant CH, spectral features described by the groups 2, 3 and 4 in paragraph 
3 . 4.1 (too narrow CH 4 bands with less good relative intensities). 

We also tested a slightly more complicated representation, R imn , by replacing the pure OL, unit 
U c of Rmcn by an intimate mixture U c+m . This did not significantly improve the fit, especially at the 

(too) strong bump around 2.28 pm. 

3.4.3 Synthesis 

The series of simulations we have performed clearly allow us to select as the best representations 
R„ and R„\. which are characterized by the presence of one or two stratified units. However, 
situations without any stratification, composed of a geographical mixture of three distinct units 
(pure CH., molecular mixture N r CH.-CO and fine grains of nitrogen), cannot be excluded even if 
they lead to less satisfactoiy fits of die observed spectrum. For both kinds of simulations we point 

out that : 

. a pure methane / diluted methane segregation appears to be necessary, whether it is vertical and 
horizontal or purely horizontal, but not granular. 
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• the areas with at least the superficial layer composed of pure methane always represent an 
important fraction of the whole surface of the visible hemisphere of Pluto ( 88 , 69 and 46 % for 
R scn , R ss ' n and Rmc n respectively). They are characterized by a granular texture of relatively small 
grains (a few hundred micrometers). 

• the concentrations of methane and carbon monoxide diluted in nitrogen are about four times 
higher on Pluto than on Triton, despite the similarities between the two bodies (Quirico et al. 
1998). 

• the huge mean grain size (of the order of several centimeters) that is obtained for the molecular 
mixture N 2 -CH 4 -CO when it forms an optically semi-infinite substratum probably is not real. 
But it indicates that these deposits present very few optical interfaces per unit volume and 
therefore probably have a compact polycrystalline texture with large evolved crystals rather than 
an accumulation of independent particles. This kind of situation was considered by Eluszkiewicz 
etal. (1998). 

V 

• A spectrally neutral bright area (or one with a very small spectral slope and/or very weak 
features) is needed over about 12% of the surface. Pure nitrogen with fine grains (order of 
10 pm) was considered here as the most straightforward candidate, but this nitrogen can in fact 
contain as much CH4 and CO as the molecular mixture unit does. Simulations with such a 
mixture for the U n unit, show that only very limited parameter changes for the 3 main surface 
representations are needed to recover very similar fits. Any other bright neutral candidate in the 
near-infrared range can be considered for this unit. 

• water ice, the most widespread compound throughout the outer Solar System, does not seem to 
be necessary to obtain very satisfactory fits of the spectrum of Pluto near the maximum of its 
visible light curve (around 200 ° longitude). 
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4. Interpretation and discussion 


The existence of seasonal transports of volatile molecules at the surface of Pluto has been 
suggested 15 years ago and first modeled by Cruikshank and Silvaggio (1980), Trafton et al. (1983) 
and (Stem et al. 1983). However at that time, only methane was identified and then considered in 
the models. The discovery of solid N 2 and CO (Owen et al. 1993) led to consideration of nitrogen 
(with a trace amount of CO) as the main volatile molecule driving the interhemispheric transport 
and methane as an « impurity » having a special behavior due to its lower volatility (Trafton 1990, 
Benchkoura 1996, Stansberry et al. 1996). As a result, a possible segregation of nitrogen and 
methane was already predicted by Lunine and Stevenson (1985) but for the surface of Triton. Using 
a globally averaged seasonal model, Trafton (1990) predicted that a layer of pure CH 4 only a few 
molecules thick should overlay the N 2 -CH 4 solid solution (« detailed balancing layer »). On the 
other hand Benchkoura (1996), using a time-dependent latitudinal model of volatile transport 
including differential sublimation of the N 2 -OLrCO mixture, obtained a geographical segregation 
of pure methane patches located at the boundary between the substrate (water ice ?) and the 
sublimating N 2 -rich polar cap. These moderately hot methane areas cover between a few percent to 
almost 40% of Pluto’s surface during the seasonal cycle. 

Finally, starting from the spectroscopic analysis of the early UKIRT/CGS4 spectra suggesting a 
possible vertical segregation of CH4 at the surface of Pluto (Schmitt et al. 1994, Quirico 1995), one 
of us performed laboratory experiments demonstrating that slow sublimation of N 2 ice with low 
CH 4 concentration (0.2%) leads to the progressive formation of a fine grained CR, layer on top of 
the sublimating N 2 -rich molecular mixture (see Stansberry et al. 1996). Extrapolating this 
experiment to Pluto’s surface, Stansberry et al. (1996) considered that the continuation of the 
segregation process may result in an ever increasing thickness of the layer of CH 4 grains, finally 
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leading lo thick lag deposits with higher temperatures than the volatile N2-rich molecular mixture. 
Modeling the vertical exchanges of methane between the surface and the atmosphere, they 
concluded that only a few percent of the surface covered by nearly pure and warm (-53 K) methane 
patches is sufficient to produce the high methane atmospheric abundance (1%) required to explain 
its elevated temperature of about 100 K (Yelle and Lunine 1989, Lellouch 1994, Strobel et al. 
1996). 

Considering the different results presented in this paper together with some considerations from 
the volatile transport models described above, we propose a ‘scenario’, describing the different 
transport and segregation phenomena that affect the identified ices and which lead to the surface 
state as revealed by the analysis of the UKIRT-CGS4 near infrared observation (see Fig. 14). 
Inversely, the following scenario coupled with some numerical results from our spectrum modeling 
will allow us to better constrain, or discard, some processes included in the different transport 
models. 

V 

As witnessed by the occurrence of methane segregation on Pluto, differential sublimation and / 
or chemical diffusion seem to affect on one side the CH4 molecule and on the other side N2 and CO 
within the volatile deposits N2-CH4-CO when subjected to solar radiation. At the surface, these 
processes create a residue of nearly pure methane forming a crust of variable thickness and divided 
into small grains by the escape of the N 2 and the CO gases or by the diffusion process. After they 
have been mixed and transported through the atmosphere, the gases settle as a molecular mixture by 
direct condensation on the coldest areas. These areas, although probably extended, cover quite a 
small fraction of the observed hemisphere, because the telescope essentially receives the light 
reflected around the subsolar point, i.e. from the well illuminated and thus sublimating surfaces. 
Since they are in thermodynamic equilibrium with an atmosphere containing the same molecules, 
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the deposits evolve more or less rapidly by metamorphism to a compact polycrystalline texture. 
Finally a new segregation cycle starts when they undergo the effects of the sun again as the seasonal 
cycle lakes place, temporally modulating the sunlight received at a given latitude. This general 
scenario can be made more specific if we look at the details of the results obtained for the different 
units of the three main representations we retained. 

First, the additional complexity of the ‘stratified representations , Rsc„ and R„~ n compared to the 
‘geographical’ one R mcn , appears to be more consistent, both at the spectroscopic and the physical 
level, with the observations and the transport models. Indeed, over most of the surface, the pure 
methane crust overlying the N 2 -CH 4 -CO mixture may be thin enough to let the radiation reach the 
compact sublimating deposits. Using the results of the model of Benchkoura (1996), we calculate 
that the sublimation of a mixture with 0.4% CH 4 leads to a methane oust about 100-200 pm thick 
over most the sublimating area. Such a layer is clearly optically thin over most of the spectrum. 
Second, the contribution of the areas displaying optically thick CH* residues, wherever they are 

f 

situated at the surface or at depth, also appears to be necessary. They attest to the existence of old 
surfaces that have sublimated for a long time (possibly during previous seasonal cycles), and 
eventually became covered over by small amounts of fresh deposits of the molecular mixture N 2 - 
CH 4 -CO (R SJ ~ n case). We note that all our models need much bigger grains for these thick crusts 
than for the methane fine crust. This may result from metamorphism processes occurring in these 
warmer CH 4 patches after the sublimation and diffusion of the underlying N 2 and CO has stopped. 
This is another proof of their age. Such a situation is particularly clear in the case of the R $s ~ n 
representation where the grain sizes differ by one order of magnitude. When it exists (R ss " n case), 
the superficial N 2 -CH 4 -CO frost layer consists of grains with sizes 20 limes smaller compared to 
the grains of the active sublimating NrCHU-CO areas. This may effectively give evidence of the 
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freshness of this thin deposit relative to the sublimating deposits, probably condensed several tens 
of years ago and strongly evolved by metamorphism and densification. The very fine grained N 2 - 
rich area should also have been deposited very recently. These grains may be linked with the 
micron sized grains that have been inferred for the thin clouds observed in the atmosphere of Triton 
(Pollack et al. 1990). 

The existence of all these differentiation and stratification processes on Pluto surely arises from 
the larger concentrations of methane and carbon monoxide that we obtain for Pluto compared to 
Triton (Quirico et al. 1998). Indeed, such phenomena clearly display a limited extent or do not exist 
at all on Triton. On Pluto, they imply the presence of an important amount of superficial pure 
methane which is certainly linked, as modeled by Stansberry et al. (1996), with the detection of 
high CH4 concentrations in the atmosphere ( ~ 1%) (Young 1994). Although the thin CH4 crust 
probably did not contribute to the upward CH4 mass flux necessary to maintain the atmospheric 
concentration in this model, there are large enough thick CH4 deposits (surface proportion - 20%) 
in both stratified representations to comply with these calculations. Indeed only 1 to 3% of surface 
cover by hot CH4 patches is enough. TTie non-detection of water ice which is only confirmed at the 
moment for the longitudes around the maximum of the light curve as well as the absence of CO 2 
ice, is another particularity of Pluto compared to Triton (Cruikshank et al. 1997, Quirico et al. 
1998). As water ice is expected to be present as a rocky lithosphere by the bulk composition models 
(see e.g. Stem 1998), the absence of such bare surface units may mean that the amount of volatile 
materials, or of their derived products, is quite large on the surface. The thick CH 4 crust may also 
hide the H 2 0 and CO 2 ices on Pluto, while the complete sublimation of the volatiles on Triton 
uncovers them. 
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Now we can determine what the best surface representations we have derived from the analysis 
of the UKIRT-CGS4 near infrared data would imply for the appearance of Pluto in the visible if it 
was imaged by the HST as performed by Stem er al. (1997). The methane, whether it is pure or 
diluted, does not present any noticeable signature below 0.62 pm, but maybe an extremely weak 
residual absorption due to overtone and combination bands. Nitrogen and carbon monoxide, which 
are quite featureless in the near infrared, are totally transparent in the visible. Consequently, a 410 
nm HST image of the Pluto’s disk as described by the R SC n, R*s n or Rmcn representations would 
show a quite bright surface on the average with some slightly darker regions corresponding to the 
compact polycrystalline N2-CH4-CO sublimating deposits. Indeed for the latter area, the photons 
travel over a much longer path length before they can exit from the surface and have then more 
chance to be absorbed by the residual absorption of CH4. The problem is that the dynamic range of 
albedo across the planet locally exceeds a factor of five ( 25 % when disk-averaged) on the real HST 
images (Stem et al. 1997). Such a contrast cannot be explained by the exclusive occurrence of 
highly translucent materials like the evolved deposits of NrC^-CO. The CCD spectra of Pluto 
recorded at various rotational phases also display a systematic decrease of albedo toward the blue 
(red slope) (Buie and Fink 1987, Fink and Di Santi 1988, Grundy and Fink 1996) that cannot be 
explained by the ice absorption so far. Nevertheless photochemical processes can affect these basic 
molecules. They may be particularly intense within the superficial crusts of pure methane, first 
forming ethane, possibly identified through the 1.75 pm band, next larger organic molecules and 
then leading to a very complex mixture with some compounds strongly absorbing in the visible. 
The results of the HST, the CCD and UKIRT-CGS 4 observations provide useful constraints on the 
optical properties of these unknown compounds. Their absorption coefficient should present quite 
high values in the visible progressively decreasing around 1.0 pm from where it should have weak 
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values in the NIR, contrasting by a factor of 100 at least with the value at 410 nm. Different kinds 
of materials may be suitable : some hydrocarbons and nitriles displaying a high carbon to hydrogen 
ratio (Gougeon 1998), like organic tholins (Khare el al. 1984) or kerogen for example. But minerals 
can also be considered. For most organic molecules, the existence in the near-infrared of many 
bands of various intensity, which are unobserved, implies that their proportions on the surface of 
Pluto must be less than a certain threshold, maybe of the order of 0.001%, given the absorption 
level of methane in the near infrared continuum and the S/N of the spectrum. More complex 
mixtures, like tholin, as proposed by Gnindy and Fink (1996), have optical constants consistent 
with our constraints. They may either replace nitrogen in the neutral unit of our representations or 
they may be mixed with methane in the old and thick CHL, crusts. This problem should be 
investigated more thoroughly by conjointly modeling the visible and the near-infrared spectra. So, 
we can note that the search for the dark components on Pluto, even if it is a bit constrained, remains 
an open question. 

V 

5. Conclusion 

The use of a polyvalent and efficient modeling tool, applied to the analysis of high quality 
spectrophotometric data, allowed us to demonstrate for the first time that on the surface of Pluto 
there exist geographical and vertical segregations of solid methane from the other detected ices, 
N 2 and CO. This implies the presence of extended surfaces with superficial pure methane where an 
intense photo-chemistry may take place, leading to the dark material evidenced by the visible HST 
images. Additionally, such methane rich surfaces give a possible explanation for the large 
abundance of CK, in the atmosphere . All the ices, N 2 , CH,, and CO seem to condense, evolve 
and sublimate in a complex manner during the seasonal cycle. The kind of horizontal and vertical 
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distributions we have obtained, at least for one hemisphere, allowed us to give a specific 
description of these processes. 

This work only represents the first part of an important series of analyses to be performed in 
the near future on observations that correspond to different rotational phases of Pluto. Such 
analyses will certainly reveal other interesting information, in particular concerning the 
geographical localization of the different surface units distinguished by the modeling of the 
infrared spectra. As a consequence, it will be easier to study the link between these units and the 
albedo features of the HST visible images whose chemical and physical nature could be then less 

ambiguously determined. 
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Tables 


Table I 

The different parameters of the basic surface representations R,, R, and R* 


First basic representation Rj 

Second basic representation R s 

Third basic representation Rg 

(granular mixture) 

(stratification) 

(geographical mosaic) 

G c [CH4] : 

U [CH4J : 

U c [CH4] : 

• number proportion p c of 

• optical depth t c 

• surface proportion c P 5urf c 

pure CH4 grains 

• mean grain size r c 

• mean grain size r c of pure 

• mean grain size r c 

• scattering parameter gc 

CH4 

• scattering parameter gc 


• scattering parameter gc 

G m [N2-CH4-CO] : 

Lm [N2-CH4-CO] : 

U m [N2-CH4-CO] : 

• number proportion p m of the 

• concentrations Cch 4 and Cco 

• surface proportion (1 -<?*"“*„,) 

molecular mixture grains 

of CH4 and CO in the N 2 

• concentrations Cch 4 and Cco 

• concentrations Cch 4 and Cco 

matrix 

of CH4 and CO in the N 2 

of CH4 and CO in the N 2 

• mean grain size r m 

matrix 

matrix 

• scattering parameter gm. 

• mean grain size r m 

• mean grain size r m 


• scattering parameter g m 

• scattering parameter g m 
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Table II 


Three unit surface representation R^n and its best fit parameter values 











Table III 


Three unit surface representation R*"„ and its best fit parameter values 



Ui [CH4] : Utf [N2-CH4-CO] : [N 2 ] : 

• optical depth Xci = 0.5 • optical depth %„& = 0.8 • mean grain size r„ = 10 pm 

• mean grain size r c i = 200 pm • concentrations • scattering parameter g„ = 0 

• scattering parameter g c i = 0 C m2 cH4 = 2.0% ; C™ co 

0 . 2 % 

• mean grain size 

rm 2 = 3.0 10 3 pm 

• scattering parameter g^ = 0 
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Ln,, [N 2 -CH 4 -CO] : 

Ui [CH 4 ] : 

• concentrations 

• mean grain size 

C m, cH 4 = 0.42% ; CT'co = 

r c2 = 1500 pm 

0 . 2 % 

• scattering parameter gc 2 = 0 

• mean grain size 


r mJ = 6.0 10 4 pm 


• scattering parameter g m i = 0 
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Table IV 


Three unit surface representation R^,, and its best fit parameter values 



[N2-CH4-COJ: [CH4] : [N 2 ]: 

• concentrations • mean grain size r c = 200 pm • mean grain size r„ = 10 pm 


Cch 4 = 0.3% ; Cco = 0.1% • scattering parameter gc = 0 • scattering parameter g„ = 0 

• mean grain size 

r m = 9.5 10 4 pm 

• scattering parameter g m = 0 
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Figure captions 


FIG. 1. : Visible geometric albedo light curve of the Pluto-Charon as calculated by Buie et al . 
(1997). The location of Pluto (sub-Earth longitude) during the three 1995 UKIRT/CGS4 
spectrophotometric observations, the 1992 and 1994 observations used to check our calibration 
operation, as well as at the time of the eclipse observation of Charon in 1987, are indicated. 

FIG. 2. : Calibration operation applied to the 05/11/95 observation of Pluto. The H region (dot- 
dash line) is adjusted to the K region by a multiplication factor to form a new coherent spectrum 
(solid line). This spectrum is compared to better calibrated spectra from 05/27/92 and 04/16/94 at 
similar rotational phases (dotted and dashed lines). 

FIG. 3.: Low resolution modeling (solid line) of the spectral geometric albedo of Charon 
calculated by Roush et al. (1996) (dashed line), from the 04/23/87 mutual event observations (Buie 
et al. 1987). In this model 60% of the surface is covered by H 2 0 ice (grain size 50 pm) and the 
remaining by a spectrally neutral component. 

FIG. 4. : Geometric albedo of the Pluto-Charon system from the 15 May observation (dotted 
line), of Charon (dashed line) at high resolution calculated by fitting the low resolution 
observations (see FIG. 3), and of Pluto alone (solid line) after extraction of the Charon 
contribution to the total flux. 

FIG. 5. : Identification of the absorption bands and illustration of the choice of the significant 
characteristic spectral features on the spectrum of Pluto that are good criteria for the estimation of 
the modeling quality, a) H region, b) K region. 
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FIG. 6 . : Comparison of the 1995 albedo spectra (a: H region; b: K region) of Pluto with laboratory 
transmission spectra of pure CPL, ice at 40 K (band positions indicated by vertical dotted lines) and 
of CH 4 diluted in P-N 2 ice at 40 K (positions : vertical dot-dash lines). The spectrum of Triton 
(09/07/95) is also shown for comparison (from Quirico et al. 1998). All spectra are shifted 

vertically for clarity. 

FIG. 7. : Two ‘best’ simulations for a simple geographical mosaic of a neutral unit U n (fine grained 
N 2 ) and a molecular N 2 -CH 4 -CO mixture (R mn ) or pure CH 4 (Ren). 

FIG. 8 . : The ‘best’ simulation for a simple geographical mosaic (Rj n )of a neutral unit U n (fine 
grained N 2 ) and an intimate (granular) mixture of pure CR, grains and N 2 -CH 4 -CO molecular 

mixture grains (U c+m ). 

FIG. 9. : Surface representation R scn : the best fit of this representation to the spectrum of Pluto. 
The best fit parameter values are reported in Table n. < 

FIG. 10. : Surface representation R ss ~„ : the best fit of this representation to the spectrum of Pluto. 
The best fit parameter values are reported in Table m. 

FIG. 11. : Surface representation R mC n : the best fit of this representation to the spectrum of Pluto. 
The best fit parameter values are reported in Table IV. 

FIG. 12. ; Comparison between the simulations corresponding to the representations Rscn and R ss n 
(see text), a) H region, b) K region. 

FIG. 13. : Comparison between the simulations corresponding to the representations R, 5 ~n and Rmcn 
(see text), a) H region, b) K region. 
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FIG. 14. : A possible 'scenario’ of the sublimation and re-condensation processes which affect the 
ices and the surface state of Pluto, as suggested by our modeling. 
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